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The Vlasov equation with strong magnetic field 
and oscillating electric field as a model for isotop 
resonant separation * 

Emmanuel Frenod & Frederique Watbled 



Abstract 

We study the qualitative behavior of solutions to the Vlasov equation 
with strong external magnetic field and oscillating electric field. This 
model is relevant to the understanding of isotop resonant separation. We 
show that the effective equation is a kinetic equation with a memory 
term. This memory term involves a pseudo-differential operator whose 
kernel is characterized by an integral equation involving Bessel functions. 
The kernel is explicitly given in some particular cases. 

1 Introduction 

This paper gives a mathematical analysis of a model related to isotop resonant 
separation. We undertake this model using homogenization methods applied to 
the Vlasov equation. 

The mathematical framework of this paper is the investigation of the in- 
fluence of oscillations generated by strong fields in the Vlasov equation. It 
completes the works led in Frenod [TT], Frenod and Hamdache [H], Frenod 
and Sonnendriicker [TH [151 [Ml HZl [H], Frenod, Raviart and Sonnendriicker 
[13], Golse and Saint Raymond [19l[20], Saint Raymond [37], Brenier [5], Gre- 
nier [22l [23], Jabin [25], Schochet [40], Joly, Metivier and Ranch [26]. We 
also refer to mathematical and physical works where similar methods are used: 
[Ml[271llHlllIllMll3ni[inill[MlllI]- The goal here is to exhibit the efl'ect of 
the interaction of the oscillations induced by a strong magnetic field with the 
oscillations of the electric field. In the case when both oscillation frequencies are 
the same, resonant phenomena appear leading to memory effects in the effective 
equation. Concerning this topic of memory effects, this work is an important 
step in the understanding of non local homogenization previously analyzed in 
Sanchez-Palencia [3S], Tartar [151133], Lions [21], Amirat, Hamdache and Ziani 
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[31 2] and Alexandre [T], in the sense that we are able, here, to exhibit explicit 
memory terms for a physically relevant problem. 

From the physical point of view, this study contributes to the understanding 
of phenomena appearing during isotop resonant separation experiments. Isotop 
resonant separation consists, in a plasma made of several kinds of ions, in heating 
only the ions of one given species. Those ions can be then easily extracted from 
the plasma by ad hoc devices. In order to reach this goal, the plasma is, inside a 
cylinder (with length about 1 m and radius about 10 cm), submitted to a strong, 
static and homogeneous magnetic field B (about 0,4 T). Under its action, each 
particle moves helicoidally around the magnetic lines with pulsation 

_q\B\ 



called cyclotron pulsation, where q and m stand for electric charge and mass of 
the considered particle. Let rh and q be the mass and the charge of the ions to be 
heated. If an electric field E, oscillating with pulsation is moreover applied 
to the plasma, the particles of the considered isotop are resonating with it and 
then acquire energy. For a detailed description of isotop resonant separation, 
we refer to Omnes |36] , Louvet and Omnes [34] , Schmitt 39J , Dawson et al. [9] , 
Compant La Fontaine and Louvet [SI [7] . 

As a model enabling us to understand some aspects of isotop resonant sep- 
aration, we introduce a small parameter e and consider the Vlasov equation 

dtf + V ■ + {Eit, ^-,x)+vx {B{t, ^, x) + ^)) ■ V„r = 0, 

f{0,x,v) = fo{x,v), 

where =^ is a strong magnetic field, M. — ei denoting the first vector of the 
canonical basis (61,62,63) of R^, and where E{t,^,x) and B{t,^,x) are fast 
oscillating electric and magnetic fields. In this equation = f^{t,x,v)] t G 
[0, r[, where T < +00, is the time, x = (xi, X2, X3) stands for the position and 
V = {vi,V2-,V3) is the velocity. This equation models the evolution of the ions 
to be heated without taking into account their interactions and the interactions 
with the other ones. This is not completely unreasonable since the ion density 
met in isotop separation experiment is relatively low (about 10^*^ ions.cm"^). 
Yet, a model taking into account self induced forces (like Vlasov-Poisson) will 
be considered in a forthcoming paper. 

We introduce the notations = Ei^ x and Q = [0,T[xn. The initial 
data satisfies 

/o>0, 0< / f^dxdv< 00. (1.2) 
Jn 

The electric and magnetic fields E{t, t, x), Bit, t, x), are C°° and 2tt periodic in 
T. Under these assumptions we have the classical a priori estimate that there 
exists a constant C independent of e such that the solution of the Vlasov 
equation (|1.1[) satisfies 

II/'^IIl~(0,T;L2(J7)) < C. 
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We deduce that up to a subsequence still denoted by e, 

f^f m L°°{0,T-,L^{9)) weak-*. 

The aim of this paper is to find out an equation satisfied by the limit /. Let us 
introduce some more notations: for any vector V in M'^ we denote by V\\ = ViCi, 
VI = V2e2 + VsCs, respectively, the parallel and perpendicular components of V 
with respect to ei. If r is in [0, 27r], we denote by r{V, r) the image of V under 
the rotation of angle r around the axis M. = ei and we define: 



E{t,x) = ^ f r{E{t,T,x),T)dT, 
B{t,x) = ^J^ r{B{t,T,x),T)dT. 



We shall denote by Jo (respectively Ji) the Bessel function of order zero (re- 
spectively of order one): 



Jo{z) — — cos(z sin r)(ir — — 

Jq Jo 



cos(— z cos r) dr, 



Ji{z) = — / cos(z sinr — r) dr = — Jq(z). 

TI" Jo 

About definitions and properties of Bessel functions, see for instance [29] or [45] , 
To finish with the notations let us precise that we shall use the notation of Taylor 
([H]) for pseudodifferential operators: 

K{z-D)h{z)^ —\— ( K{z,ik)Th{k)e''-^dk. 



(27r)3/2 

Now we are ready to state our main result: 

Theorem 1.1 (Main result) We assume that fo{x,v) depends only on x, vp 
and the modulus ofv±, and satisfies hi .2(1 . We assume also that E{t,x) does 
not depend on x\\, that — and that B ^ 0. Then the sequence (/^) of 
solutions of lil.l]) satisfies, for any T > 0, 

^ / m L°°(0, T; L^{VL)) weak-*, 

and f formally satisfies the following partial integro differential equation 

dtf{t,x,v)+vi-dxj{t,x,v)^ / K{s,t,y;D)f{s,x,v) ds 



(1.3) 

/(0,a;,w) = fo{x,v), 

where K{s,t,y; D) is a pseudo- differential operator parametrized by {s,t) in 
[0, and y — {vi,X2,X3) in MP, acting on functions of z = {xi,V2,V3) vari- 
ables. The kernel 

K{s,t,y;ik) = exp{-ikiyi{t - s))K{s,t,y;ik) (1.4) 
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where K{s,t,y]ik) is the unique solution of the Volterra equation 



Jo{\kj_ I \L{s, a, y)\)K{a, t, y, ik) da 



L{s,t,y)-E{t,y) 
\L{s,t,y)\ 



\k^\Ji{\k±\\L{s,t,y)\) (1.5) 



with 



L{s,t,y) 



E{a, y)da 



E±{a,y)da. 



(1.6) 



Remark 1.1 The assumption on variable-dependancy as the one on E\\ we 
make in the heading of the Theorem is to ensure a transversahty property re- 
quired in order to apply non local homogenization methods. Nevertheless, from 
physical point of view, as original repartition functions are often maxwellian, 
considering that /o depends only on x, and the modulus of v± is relevant. 
The two facts £'|| = and B = can be easily realized experimentally, of course 
as soon as the self induced fields are neglected. 

Remark 1.2 Equation (|1.3p does not contain oscillations. It only contains the 
mean effect, via the memory term, of the oscillations contained in equation 
(|l.ip . According to our knowledge, it is the first time that such a mean model 
is exhibited to describe the dynamic of ions during isotop resonant separation. 

In some cases, we may obtain explicit expressions of the kernel K: 

Theorem 1.2 In the particular case where 

Ej_{t,T,xj_) ^ £{t,x±)g{T,x±), 

with £ ]S.-valued and under the assumptions of Theorem ll.ll the kernel is exactly 
equal to 

K{s,t,y,ik) 



Remark 1.3 As it is stated in Alexandre [Tj, it is immediate to check that the 
kernel K appearing in p.7p satisfies: 



exp{-ikiyi{t - s)) £{t,y)£{s,y) \k±\\g{y) 



Ji{\k^\\9{v)\rs£{<^,y)d<j) 

Sl£{(J,y) da- 



(1.7) 



where 
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The following Theorem is a simplified version of the previous one, but as its 
proof is much simpler; we state it separately. 

Theorem 1.3 In the very particular case of an electric field independent of 
time: 

E^ = E^{t,x^) (1.8) 
and under the assumptions of Theorem \l.l[ we have: 

K{s, t, y, ik) = - eM-^km{t - s)) '^'^'^''"^'-^^^'^^ ~ ■ (1.9) 

t — s 

The proof uses the notion of two scale convergence introduced by N'Guetseng 
[35] and Allaire [2]. Their result is the following: 

Theorem 1.4 (N'Guetseng and Allaire) // a sequence [f^) is hounded in 
L°°(0,T; L'^{n)), then there exists a 2Tr-periodic in r profile F(t,T,x,v) in 
L°°(0,T; L°°{S.r; L'^{fl))) such that, for every 'ip{t,T, x,v) regular, compactly 
supported with respect to {t, x, v), and 2Tr-periodic with respect to t, we have, up 
to a subsequence, 

f^ip^ dt dx dv ^ / Fip dr dt dx dv, 
Q Jq Jo 

where ip^{t, x, v) = 'ip{t, |, w). The profile F is called the 27r -periodic two scale 
limit of and the link between F and the weak-k limit f is given by 

F{t, T, X, v) dr = f{t, X, v). 

It has been used by Frenod and Sonnendriicker [15l [TTl [18] in the context of 
homogenization of the Vlasov equation but under an assumption of strong con- 
vergence of the electric field. Using the same ideas we obtain first the following 
result concerning the two-scale limit of : 

Theorem 1.5 The sequence (/^) of solutions of two scale converges to- 

wards the 2TT-periodic in t profile F which is the unique solution of 

drF+{v X M) ■ VyF = 0, 
dtF + vi ■ d,,F + {r{E, -t) + vx r{B, -r)) • = ^q) 

F{Q,T,x,v) = —fo{x,r{v,T)). 

Then following Alexandre [1] we apply a Fourier transform J-' to obtain an 
ordinary differential equation satisfied by TF. Using the fact that 



^27r 

Tf{t,x,v) = 



/ !FF[t,T,x,v) dr, 
Jq 



we obtain an ordinary differential equation satisfied by J^f. Eventually we 
introduce a Volterra equation as in Tartar [42l [43] and Alexandre [1] to obtain 
our result. 



6 



The Vlasov equation with strong magnetic Geld 



EJDE-2002/?? 



2 Scaling and qualitative study 

We exhibit the important parameters playing a role when charged particles are 
submitted to a strong magnetic field and in view of which we provide the scaling 
leading to equation (jl.ip . 

Before any scaling the evolution of the repartition function f{t,x,v) rep- 
resenting at each time t the particle density standing in x and moving with 
velocity v, is given by the Vlasov equation 

dtf + v-V^f+^{E + vxB)-V^f^O. (2.1) 
m 

We define now characteristic scales: t is the characteristic time, L the character- 
istic length and v the characteristic velocity; and rescaled variables: t' , x' and v' 
hy t — It' , X = Lx' , V — vv' . We also define scaling factors for the fields E and 
^ and rescaled fields E'{t\x') and B'{t',x') by EE'{t',x'_) = E{tt',Lx') and 
BB'{t\x') = B{tt' , Lx'). Lastly, defining a scaling factor / for the repartition 
function, without forgetting that / is a density on the phase space, we define 
the rescaled repartition function /' by 

ff'{t',x',v') = L^lff{tt',Lx',vv'). 
The new repartition function is the solution to 

dt-r + ^v' ■ V,,/' + {^E'if, x') + '^v' X B'{t', x')) . V.,/ = 0. (2.2) 

Let us introduce two parameters having an important physical signification: 
uJ^ = ^ is the characteristic cyclotron pulsation and ol = = the characteristic 
Larmor radius. Looking at (12. 2p in view of these parameters, we get 

dt'f + t—,^v' ■ V,'/' + {t—,^E'{t', x') + t—,v' X B'it', x')) ■ V„./' ^ 0. (2.3) 
L vB 

Now, introducing the small parameter e, we set 

^ = e and uJ^t — -, (2.4) 
L £ 

which means that the observation length scale and the observation time scale are 
respectively large in front of the Larmor radius and the cyclotronic period. This 
regime is relevant to describe the global behaviour of the considered particles. 
We also assume that the electric force is much smaller than the magnetic one. 
It reads: 

^-e. (2.5) 
vB 

Now, in order to model the fact that the particles are submitted to an oscillat- 
ing electric field, we set that E' writes E'(t', j,x'). Concerning the magnetic 
field we assume that it is made of a constant field perturbed by an oscillating 
one, it gives B' = M + eB"{t', ^,x'). Using those assumptions and removing 
subscripts " and ', equation (|2.3p leads to (|l.ip . 
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Remark 2.1 The goal of this remark is to show that the considered scaling 
corresponds to physical situations. If the ions are potassium and the magnetic 
field magnitude is IT, the characteristic cyclotron pulsation is about 10^ (= ZI^). 
In experimental situations it is realistic to consider that the ions stay 10~^ s in 
the device (we make this value as reference time t). Beyond this, considering 
that the thermic velocity (we choose as characteristic velocity v) of the ions 
is something like lO^m.s"^, the size order of their Larmor radius is the mm 
(ol 10~^m). As the device is about one meter long, we have 

^-10-3, zu2l^lO\ (2.6) 
L 

Yet if the magnitude of the electric field is about IKm^^ we also have 

^ = 10-'. (2.7) 

As stated in Theorem ll.il the effective behaviour of equation (|1.3p involves 
memory effect. In order to give a way to understand why, we propose to begin 
by studying the following simple problem: analysing the movement of a particle 
(with mass and electric charge equal to 1) with no velocity component parallel 
to M. — ei, and submitted to the magnetic field M je. Under the mere action of 
this magnetic field, the particle rotates around M with pulsation 1/e. In other 
words, its velocity and position write 

V{t) = r(z;o, ~ J), X(t) = + e{r{va, f - |) + t;o x M), 

where vq and xq are velocity and position of the particle at t (we assume 
vq-^A — 0). As £ goes to 0, X{t) tends to xq. Yet, V(i) drives the circle \v\ = \vq\ 
faster and faster. Consequently we could say that the 'particle occupies the whole 
circle \v\ = \vo\ and has forgotten its initial direction If now an oscillating 
electric field, writing for instance 



E'it)^rie2,~*-)= \ cos(t/e) 





sin(</£) 



is applied to the particle in addition to the magnetic field, the movement of 
the particle is no more a rotation, but a spiral around the magnetic field M. 
Indeed, the considered electric field does not modify the angular velocity of the 
particle but only the modulus |f | = + v^, since we have 

V{t)^r{vo,--)+tE'{t). 
e 

Observe that 

±\V{t)\^ = 2ivo-e2 + t). 
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If uq • 62 > the modulus \v\ of the velocity increases during each rotation; 
if Wo • 62 < the modulus decreases during the first rotations (until the time 
t = —vq ■ 62) and then it increases. Consequently, we could say here that the 
dynamic of the particle strongly depends on the initial value vq . Of course this 
dependance is kept when e goes to 0. But as we just saw the particle forgets 
its initial value as s goes to 0. Because of this contradiction we need to keep 
additional information, which is contained in the memory term taking place in 
equation (|1.3|) . 

To finish this qualitative study, we shall explain the result contained in 
Theorem 11.51 using formal asymptotic expansion. If we assume the following 
ansatz of f^ 

f{t,x,v)^Fo{t,-,x,v)+eFi{t,-,x,v) + ..., (2.8) 
e £ 

if we insert this in (jl.ip and identify the terms at each order we get, at order 0: 

drFo + {vxM)- V^Fo = 0, (2.9) 

and 

drFi + {vxM)- VyFi 

= - {dtFo + V ■ V,Fo + {E{t, T, x)+vx B{t, r, x)) ■ V,Fo) . (2.10) 

We then see that the first term Fq is nothing but the two scale limit F and that 
equation (|2.9p is nothing but (|1.10b .). The second equation (|1.10b ) is given as 
a compatibility condition on Fq — F va. order that (|2.10p has solutions. 



3 Equation satisfied by the 2-scale limit. 

In this section we follow the procedure of Frenod and Sonnendriicker [15j to 
prove Theorem 11.51 We only sketch the proof and refer the reader to [THj for 
details. 

Let {f^) be a sequence of solutions of p.ip . As the sequence is bounded 
in L°°(0, T; i^(ri)), by Theorem 11.41 there exists a 27r-periodic in r profile 
F{t,T,x,v) in L°°{0,T]L°°{Rr]L'^{i^))) such that for every ip{t,T,x,v) regu- 
lar, compactly supported with respect to {t, x, v) and 27r-periodic with respect 
to T, we have, up to a subsequence, 

/ fiP'^dtdxdv^ I [ FiPdrdtdxdv (3.1) 
Jq JqJo 

The proof of Theorem 11.51 is led in three steps: 

Step 1: First we use a weak formulation of (|l.ip in 'D'{Q) with functions 
•ip'^(t,x,v) — iplt, -,x,v) where "0 is regular with compact support in {t,x,v) 
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and 27r-periodic in r, which writes: 

/ ridti!" + v- V^i:^ + {E^ {B^ + — )) ■ Vyil>^) dtdxdv 

Jq £ 

= - fo{x,v)ip%0,x,v)dxdv. (3.2) 
Jn 

Notice that 9t (?/'") = {dtip + -drip)". Multiply ([X^ by e, then let e tend to 

and apply the two scale convergence (|3.ip to deduce that F belongs to the 
kernel of the singular perturbation appearing in (jl.ip . in other words: 

drF +{vxM)- VyF = in 2?'(M^ x M^) (3.3) 

for almost every {x,t) in [0,T[xR^. 

Step 2: Next we use the following lemma (Lemma 2.3 of [E]), which character- 
izes the kernel of the singular perturbation. 

Lemma 3.1 A junction F{t,v) G L°° (Rr , L'^ (R.^)) 2n-periodic in r satisfies 

drF + {vx M)- VvF = m V'{M.r x M^) 

if and only if there exists a function G G L^(M^) such that F{t, v) = G{r{v, r)). 

According to this lemma there exists a function G in i°°(0, T; L'^iR^ x M^)) 
such that 

F{t,T,x,v)^G{t,x,r{v,T)). (3.4) 

Step 3: We denote here rj' = K3xM3^ Q' = [0,T[xrj'. The goal of this step is to 
project equation (jl.ip on the orthogonal of the kernel we identified above and to 
pass again to the limit. In order to achieve this we build test functions belonging 
to the kernel in the following way. For every regular compactly supported 
^p{t, X, u) we consider the 27r-periodic in t fimction iplt, t, x, v) = ip{t, x, r{v, t)), 
which in view of Lemma |3 . II satisfies 

driJ + {vx M)- Vyij = 0. 

We take ip^{t, x, v) = ilj{t, ^, x, v) in the weak formulation of the Vlasov equation 
(|3.2p obtained in the first step so that 

/ /'((9t^)' + {v ■ V.^A)" + {{E + vxB)- S/y^Y) dtdxdv 
JQ. 

/o(a;, v)ip{Q, 0, x, v) dx dv. 

We let £ tend to 0, use the two scale convergence p.ip and the equality (|3.4|) to 
deduce that 

2-77 

Git, X, r{v, r)) (^dt(p + v ■ + r[E + v x B,t) ■ V„(/?) dr dt dx dv 

fo{x,v)(p{0, x,v) dx dv. 
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Now we make the change of variables u — r{v,T) and perform the integration 
with respect to r over [0, 27r]. We obtain in this way the equation satisfied by 
G: 

Lemma 3.2 The function G{t, x, u) linked to the 2tt -periodic profile F by 
is the unique solution of 

dtG + ui-d^,G+iE + ux B)- V„G = 

1 (3.5) 
G{0,x,u) = —foix,u). 

ZTT 

The uniqueness of the solution of (|3.5p enables us to deduce that the whole 
sequence f^ two-scale converges to F and, because of the link between F and 
/, weak-* converges to / (recall that J^^ F{t, r, x, v) dr = f{t, x, v)). 

To prove Theorem 1 1.5 1 we rewrite the equation (|3.5p in terms of F using the 
equality F{t, r, x, v) = G{t, x, r{v, r)) and thus obtain the equation satisfied by 
F since 

r{VyF{t, T, X, v),t) ^ VuG{t, X, r{v, t)) 

and 

{r{v, t)xB) ■ r(V„i^, r) = (w x r(B, -r)) ■ V^F. 
The Theorem is then proved. □ 

Remark 3.1 At this stage we can integrate in r over [0, 2tt] the equation 

dtF + vi ■ d^,F + {r{E, -t)+vx r{B, -r)) • V„F = 

F{0,T,x,v) = -^fo{x,r{v,T)) 



and use the equality /(i, x, v) = J^^^ F{t, t, v) dr to get the following equation 
which is satisfied by /: 

dtf + vi-d^J+{E^\+vxB^\)\/„f 

+ {E±+vx B±) ■ V,„a +{ExM + vx{Bx M)) ■ V„6 = 0, 

where a{t,x,v) = J^^ F{t, t,x,v) cos t dr, b{t, x, v) — J^^ F{t, t, x, v) sin r dr. 

4 Equation satisfied by the weak-* limit. 

In this section we prove our main results Theorems II. 1) 11.21 and 11.31 by using 
Fourier transform. Recall that under the assumptions of Thcorcm l 1 . 1 1 the 2-scale 
limit F satisfies the equation 

dtF + «i • d.,,F + r{E, -r) • V„^^ - 

1 (4.1) 

F{0,T,x,v) = — /o(a;;W||, 
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with 



r{E{t,x±),-T) 







E2{t, x±) COS T + E^{t, x±_) sinr 
— i?2 , a;^ ) sin t + E^ {t,x±) cos r 

Thanks to the hypothesis the only derivatives of F involved in the equation 
are with respect to xi, V2, V3 and the coefficients only depend on vi, X2, X3. 
The transversality assumption required for non local homogenization methods 
is realized. For convenience we rename the variables by 



and define 



Xi 

V3 



X3 



H{t,T,z,y) = F{t,T,x,v) 
h{t,z,y) = f{t,x,v), 



ho{z, y) = h{0, z, y) = /(O, x, v) = /o(a;, v), 
a{t,T, y) 



Vl 

E2 cos r + i?3 sin t 
—E2 sin T + E^ cos r 



so that (|4.ip becomes 



dtH{t, T, z, y) + ait, r, y) ■ \JzH{t, r, z, y) = 0, 
H{0,T,z,y) = -^ho{z,y). 



(4.2) 



Applying a Fourier transform in the z variable we get the ordinary differential 
equation 

dtJ-H{t, T, k, y) + ik ■ a{t, r, y)TH{t, t, k, y) — 0, 

1 (4.3) 
TH{0,T,k,y) = —Tho{k,y), 



2tt 



which has the explicit solution 



1 /■* 

TH{t, T, k, y) = —Tho{k, y) exp{-ik- / a(cr, t, y) da) 
^'^ Jo 



271 



Tho{k,y)exp{-ikiyit)exY>{-ik- / r[E{a,y), ~t) da) 



As 



we know that 



f(t,x,v) 



F{t, T, x, v) dr 



/•2tt 

Th{t,k,y)= TH{t,T,k,y)dT, 
Jo 
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hence 

1 /■^'^ /■* 

J^h{t, k, y) = J^ho{k, y) exp{-ikiyit)— / cxp [-ik- / r{E{a, y), -T)da) dr. 

Jo Jo 

We introduce the quantities 

L{s,t,y) = / E{a,y)da, (4.4) 

J S 

A{s,t,y,Z)^^ [ exp{^Z - [ r{E{a,y),-T)da)dT, (4.5) 

^'^ Jo Js 

SO that 

Th{t, fc, y) = Thn{k, y) exp(-jfcij/ii)A(0, t, y, ifc), 
and differentiating with respect to the following equation satisfied by Th: 

dtTh{t, k, y) + ikiyiTh{t, k, y) = ThQ{k, y) exp{-ikiyit)dtA{Q, t, y, ik), 
Th{Q,k,y) = Tho{k,y). 

Now we use the method of Tartar [42l [43j (see also Alexandre [I]) to form an 
integro-differential operator. We define 

C(s, t, y, Z) = ~dtA{s, t, y- Z^, Z^) ex^{-Z^y^{t ~ s)) for Z e 

so that 

dtTh{t, k, y) + ik,y,n{t, k, y) = ~J^h„{k, y)C(0, t, y, ik). (4.7) 
We denote by D{s, t, y, Z) the solution of the Volterra- Green equation 

D{s, t, V^Z)- C{s, a, y, Z)D{a, t, y, Z) da = C(s, t, y, Z). (4.8) 

We replace C(0, i, y, ik) in equation (|4.7[) by its integral form, which gives 
dtTh{t, k, y) + ikiyiTh(t, k, y) 

= — J^ho{k,y)D{0,t,y,ik) + / Tho{k,y)C{0,a,y,ik)D{a,t,y,ik)da 

Jo 

= - Tho{k,y)D{Q,t,y,ik) ~ / dtJ^h{a, k,y)D{a,t,y,ik) da 

Jo 

— ikiyi / !Fh{a, k,y)D{a,t,y,ik) da. 



Integrating by parts in the second term, observe that 

A{t, t, y, Z) = 1, D{t, t, y, Z) = -^dtA{t, t, y, Z) = 0, 
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and we obtain eventually that 

dtJ-h{t, k,y) + ikiyiJ^h{t, k,y) ~ / K{a,t,y,ik)Th{a, k,y) da, 

Jq 

Th{Q,k,y) = Tho{k,y), 

where K{s,t, y, Z) = dsD{s,t,y, Z) - ikiyiD{s, t,y, Z). 
To simplify the expression of K we define D by 

D{s,t,y,Z) = eyiY>{Ziyi{t - s))D{s,t,y,Z). 

Then replacing D in the expression of K we obtain that 

K{s, t, y, Z) = exp ( - Ziyi{t - s))K{s, t, y, Z), 

which is equation (jl.4p . with 

Kis,t,y,Z)=d,Dis,t,y,Z), 



(4.9) 



(4.10) 



(4.11) 



and replacing D in the Volterra equation 
of the Volterra equation 



we 



obtain that D is the solution 



D{s,t,y,Z) + 



I dtA{s,a,y,Z)D{<j,t,y,Z)d<j^-dtA{s,t,y,Z). (4.12) 

We observe that D does not depend on Zi . Performing an integration by parts 
we get the following equation satisfied by K = dgD: 

[ Ais,a,y,Z)K{a,t,y,Z)da^dtA{s,t,y,Z). (4.13) 

J s 

Let us now simplify the (|4.5p expression of A: for fixed s,t,y,k, we choose 
a and /3 in [0, 2tt] such that 

r{L{s,t,y),a) = \L{s,t,y)\e2, 
k±_ = |fcj_|r(e2,/3). 



Then 



1 f^^ 
A{s,t,y,ik) = — exp 
27r Jo 



TT / *3xp 



27r 
1 

1 

1 

271 

1 

2^ 



r2w 



r2w 



TT / exp 



r2TT 



TT / CXp 



COS 



[ — ik ■ j r(^E{a, y), — r) rfcr) dr 
{-ik- r{L{s,t,y), -r)) dr 
( - i\kj_\\L{s,t,y)\r{e2, P + a + t) ■ 62) dr 
( — i\k±\\L{s,t,y)\ cost) dr 
(|fc_L||L(s,t,2/)| cost) dr, 
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so that 

A{s,t,y,tk) = Jo{\k^\\L{s,t,y)\). (4.14) 
Differentiating with respect to t, we obtain 

dtA{s,t,y,tk) = - ^^'f;y^ '^^^^y^ \k^\J,{\k^\\L{s,t,y)\). (4.15) 

l-i^Ks, t, y)\ 

Hence K{s,t,y,ik) is the unique solution of the equation 



Jo{\k±\\L{s,a, y)\)K{a,t,y,ik) da 

= -%f^^T^|fcJJi(|fcJ|i(^,t,y)l) (4.16) 
\L{s,t,y)\ 

for every s,t e [0, T], that is (fT31) . 

Applying formally the inverse Fourier transform to (14. 9|) we find that 



dth{t,z,y)+yi-dz^h{t,z,y)= / K{a,t,y; D)h{a, z,y) da. 

Jo 

Now 

HtyZ^y)= H{t,T,z,y)dT F{t,T,x,v)dT ^ f{t,x,v), 
Jo Jo 

from which we deduce p.3p satisfied by /, and completes the proof of Theo- 
rem O □ 

We shall now treat the particular case of Theorem I 1 . 31 where E is independent 
of time: in this case E — E± is independant of time too and we have 

L{s,t,y) = it-s)E^{y), 

so that using the parity of Jo and the imparity of Ji the last equation becomes 



Jo{\k^\\E^{y)\{a - s))K{a,t,y,ik) da 

= ~\E^{y)\\k^\J,{\k^\\E{y)\{t ~ s)) for every s,t e [0,T]. 
For sake of clarity we fix for the moment y and k and we set 

C ^ \k^\\E^{y)\. 

It is easy to see that for every h <e]—T, +T[, the function (s, t) K{s + h, t + h) 
is solution of the equation on [—h,T~h] x [—h,T — h], so that K(s + h,t + h) = 
K{s, t) for every s,t G [0, T] n [— /i, T — h\ by uniqueness of the solution. So we 
can set 

G{t - s) = k{s,t). 
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Replacing in the equation, performing a change of variable and setting x — t — s, 
one gets 

Jq{C{x - u))G{u) du = ~CJi{Cx) for every x e [-T,T]. 

As G is easily seen to be even this is equivalent to 

Jo(C(a; - u))G(u) du = -GJi{Gx) for every x G [0, T]. (4.17) 

Performing a Laplace transformation in x gives the equality of functions 

C,h{C-) xCG = C{-CJi{C-)), 

from which we deduce the Laplace transform CG of G, and by the inverse 
Laplace transform (see for instance the formulas for Laplace transformations in 
[46] ) we obtain that 

GMGx) 
G(x) = on 0, i , 

X 

hence on [~T, T] by parity, and eventually one gets 

^ = ^ CJiiGit- s)) ^^^^^ ^^^^ ^^^^g^ 
t — s 

giving Theorem [131 D 
Concerning the growth in terms of power of fc, we see that 

\Kis,t,y,ik)\ < |fc±Pp||oo||J{||oo < (l + \k\'')\\E\U 

Now we shall prove Theorem 1 1.2 1 Here E can be written in the form 

E±{t,T,x±) = £{t,x±)g{T,x±), 

with £ M-valucd and g R^-valued. For simplicity we fix x and we write 

E^{t,T,x^) = E^it,T) = £{t)giT). 

Then E{t) — S{t)V, where V is the vector V = ^ J^^ r{^g{T),T) dr, and 



Lj^(s, t)^V £{a) da = V{^{t) - $(s)) 

J s 

where we note $ a primitive of £. The equation (|4.16p becomes 



Jo{\k^\\V\\<S>{a) - <i>{s)\)K{a,t,ik) da 
$(i) - $(s) 



|$(t)-$(s)| 



k^\\V\£{t)J^{\k^\\V\mt)-^s)\), 
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in other words 
rt 



Jo{\k^\\V\{<i>{a)-^s)))Kia,t,ik)da 

= -\k^\\vmj,{\k^\\vmt) - $(s))) (4.19) 

for every s,t [0,T]- Now if £{t) is strictly positive (respectively strictly 
negative) for every t E]a,b[, we have ^'{t) — £{t) > (respectively < 0) hence 
$ is bijective from [a, b] onto the interval / of extremities $(a), $(6). We make 
the change of variables u = ^{a) in the integral and we obtain 

M\k^\\V\{u-a))K{^-\u),^-\f3),^k) J^^^^ 

- -\k^\\V\£{'^-\P)).h{\ki^\\V\{p - a)) for every a,/? e /, 
where we have set a = <i>(s) and /3 ~ $(i)- We set 

A-(<i>"l(u), $-!(/?)) 

^"'''^ £($-i(a))f(<i>-i(/3)) 

and then we have 

Ja{\kj_\\V\{u - a))r{u, P) du 

= -|fc_L||V^|Ji(|fc±||V^|(/3-a)) for every a,f3€l. 
We met the same equation in the preceding case. The solution is 

l3 — a 

which gives 

cu.c^ ^ \k^\\v\M\k^\\v\m) ~ m)) , ,1 

K{s,t) = -£{t)£{s) ^^^^ _ ^^^^ for every s,te [a,b]. 

(4.20) 

At this point we know K{s, t) for s, t both in an interval where £ is strictly 
positive or strictly negative. For such s, t we have, replacing in (|4.19|) . the 
equality 

jQ{\k±\\V\{<^{a) - <^{s)))£{a) $(^) _ $(^) '^'^ 

= J,{\k^\\V\im^m))- 

We can write 

Jo{\k^\\V\ma) - ^s))) 

n— 
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and write the product of the two series as 

.h{\k^\\v\im-H'^))) 



Jo{\k^\\Vma)~^s)))- 



- $(cr) 

oo 

= ^P„(s,t)(<i>W-$(a)r 



Integrating term by term we get 

so that the foUowing equaUty holds: 

n+1 ■'^ n! 

n— n— 1 

As this is an equaHty of series which has nothing to do with the sign of 8 we 
deduce that the formula (|4.20p giving K{s, t) is valid for every s, t in [0, T]. The 
inequality 

\K{s,t,y,ik)\ < ||£||oo|5(y)l'(l + |fc|') 



is immediate. This completes the proof of Theorem 11.2 
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